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■ We study Gamow- Teller strength distributions of '^^Ge and '^^Se within a deformed QRPA formal- 
ism, which includes residual spin-isospin forces in the particle-hole and particle-particle channels. 
We consider two different methods to construct the quasiparticle basis, a selfconsistent approach 

' based on a deformed Hartree-Fock calculation with density-dependent Skyrme forces and a more 

^SJ I phenomenological approach based on a deformed Woods-Saxon potential. Both methods contain 

pairing correlations in the BCS approach. We discuss the sensitivity of Gamow- Teller strength 
distributions to the deformed mean field and residual interactions. 
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O ! I. INTRODUCTION 

m 
Q 

The nuclear double /3-decay process is widely considered [1] as one of the most important sources of information 
about fundamental issues, such as lepton number nonconscrvation and massive neutrinos, that can be used to test 
' . the Standard Model. 

Q Theoretically, a condition to obtain reliable estimates for the limits of the double /3-decay half-lives is that the 

^ . nuclear structure involved in the process through the nuclear matrix elements can be calculated correctly. The 
S ■ proton-neutron quasiparticle random phase approximation (pnQRPA or QRPA in short) is one of the most reliable 
^ , and extended microscopic approximations for calculating the correlated wave functions involved in (3 and double /3 

JJ"! ' decay processes. The method was first studied in ref. [2] to describe the (3 strength. It was developed on spherical 

rS \ single-particle wave functions and energies with pairing and residual interactions. 

■ The QRPA method was also successfully applied to the description of double /3-decay [3] after the inclusion of a 
" " ' particle-particle [pp) residual interaction, in addition to the particle-hole (p/i) usual channel. Many more extensions 

of the QRPA method have been proposed in the literature, see ref. [4] and references therein. 

An extension of the pnQRPA method to deal with deformed nuclei was done in ref. [5] , where a Nilsson potential was 
used to generate single particle orbitals. Subsequent extensions including Woods-Saxon type potentials [6], residual 
interactions in the particle-particle channel [7], selfconsistent deformed Hartree-Fock mean fields with consistent 
residual interactions [8] and selfconsistent approaches in spherical neutron-rich nuclei [9], can also be found in the 
literature. Nevertheless, the effect of deformation on the double /3-decay processes has not been sufficiently studied 
[10,11]. 

In ref. [8] , ground state and /3-decay properties of exotic nuclei were studied on the basis of a deformed selfconsistent 
HF-I-BCS-I-QRPA calculation with density dependent effective interactions of Skyrme type. This is a well founded 
approach that has been very successful in the description of spherical and deformed nuclei within the valley of stability 
[12]. In this work we extend those calculations to the study of the dependence on deformation of the single (3 branches 
that build up the double /3 process. We focus on the example of the double /3-decay of ^^Ge and study /3^ Gamow- 
Teller (GT) transitions to the intermediate nucleus as well as the /3+ Gamow- Teller transitions of the daughter nucleus 
^^Se to the same intermediate nucleus. We discuss the similarities and differences of using different single particle 
mean fields of Woods-Saxon (WS) and Hartree-Fock (HF) types. 
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In sect. II, wc present a brief summary containing the basic points in our theoretical description. Section III contains 
the results obtained for the bulk properties of ^^Ge and ^®Se and a comparison of our results to the experimental 
available information. In sect. IV we present our results for the GT strength distributions and discuss their dependence 
on the deformed mean field and residual interactions. The conclusions are given in sect. V. 

II. THEORETICAL APPROACH 

In this Section we describe the QRPA formalism used in this work, which is based on two different assumptions 
for the deformed mean field, a Woods-Saxon potential and a selfconsistent mean field obtained from a Hartree-Fock 
procedure with Skyrme forces. 

In the first approach we use a deformed WS potential with axial symmetry to generate single particle energies and 
wave functions. The parameters of this potential are taken from the work of Tanaka et al. [13]. This parametriza- 
tion was proposed originally for spherical nuclei ranging from ^^O to ^°^Pb, but the derived isospin dependence of 
the parameters allows an extension to deformed nuclei as well. Previous QRPA calculations have shown that this 
parametrization provides realistic levels also for deformed nuclei and good results on Ml excitations were obtained 
[14] for nuclei in various mass regions as well. 

In these calculations, the quadrupole deformation of the WS potential 02 is usually determined by fitting the micro- 
scopically calculated ground state quadrupole moment to the corresponding experimental value. The hexadecapole 
deformation /3i is expected to be small for these nuclei and wc assume it is cqiial to zero. 

On the other hand, we also perform selfconsistent microscopic calculations based on a deformed HF method with 
density-dependent Skyrme interactions. We consider in this paper the force Sk3 [15] and the force SG2 [16] that 
has been successfully tested against spin and isospin excitations in spherical [16] and deformed nuclei [8,17]. For the 
solution of the HF equations we follow the McMaster procedure that is based in the formalism developed in ref. [18] 
as described in ref. [19]. Time reversal and axial symmetry are also assumed here. 

In both schemes, WS and HF, the single-particle wave functions are expanded in terms of the eigenstates of an 
axially symmetric harmonic oscillator in cylindrical coordinates, which are written in terms of Laguerre and Hermite 
polynomials. The single-particle states \i) and their time reversed \i) are characterized by the eigenvalues fl oi Jz, 
parity tt,, and energy Cj 

l^)=E 2 ' S 17Vn.n,AI]) (2.1) 

with Oj = A -I- S > |, and 

N)=E^^^4^ E C'k.„.AE(-l)^"''l^nrn.-A-S) (2.2) 

JV nr,n^,A>0,E 

with fli = —Qi = — A — E < — |. For each N the sum over nr,nz,A > is extended to the quantum numbers 
satisfying 2nr + Uz + A = N. The sum over N goes from = to = 10 in our calculations. 

Pairing correlations between like nucleons are included in both cases in the BCS approximation with fixed gap 
parameters for protons A^r, and neutrons A^. 

The number equation in the neutron sector reads 

2E^' = ^ (2.3) 

i 

where vf are the occupation probabihties 

; = (2-4) 

in terms of the quasiparticle energies 

Ei = ^(ei-A.)'+A2 (2.5) 

These equations are solved iteratively for the WS and HF single-particle energies to determine the Fermi level A,^ and 
the occupation probabilities. Similar equations are used to determine the Fermi level and occupation probabilities for 
protons by changing N into Z, Ay into A^r, and into A,r- 



A. 
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The fixed gap parameters are determined phenomenologically from the odd-even mass differences through a sym- 
metric five term formula involving the experimental binding energies [20] : 

A, = l[B{N -2,Z)-4B{N-1,Z) + 6B {N, Z) 
o 

-AB{N + 1,Z) + B{N + 2,Z)] (2.6) 

A similar expression is found for the proton gap by changing by Z and vice versa. For ''^^Ge we obtain = 1.54 
MeV, = 1.56 MeV and for ^^Se we obtain = 1.71 MeV and A^ = 1.75 MeV. 

Therefore, at the quasiparticle mean field level, we can observe several differences with respect to the treatment 
of the mean field in terms of HF or WS potentials. The most important is that the quadrupole deformation of the 
ground state is determined selfconsistently in HF and no explicit input parameter is needed. Other differences come 
from the structure of the two-body density-dependent Skyrme force that contains terms absent in the WS potential, 
such as a spin-spin interaction in the selfconsistent mean field through the spin exchange operators of the Skyrme 
force. 

Now, we add to the mean field a spin-isospin residual interaction, which is expected to be the most important 
residual interaction to describe GT transitions. This interaction contains two parts. A particle-hole part, which is 
responsible for the position and structure of the GT resonance [7,8] and a particle-particle part, which is a neutron- 
proton pairing force in the J'^ = 1+ coupling channel. 

V§^ = 2x'g% E /3+=^(Ha;,|7r)a^a.; (2.7) 



z.)a+a±. (2.8) 

K Trv 

The two forces ph and pp are defined with a positive and a negative sign, respectively, according to their repulsive 
and attractive character, so that the coupling strengths x ^nd k take positive values. 

The particle-hole residual interaction could in principle be obtained consistently from the same Skyrme force used 
to create the mean field as it was done in refs. [8] to study exotic nuclei. However, in this paper we use as a first 
attempt the coupling strengths from ref. [7]. In this reference, the strengths XxjT^ Kq'ji are considered to be 
smooth functions of the mass number A, proportional to . The strength of the ph force is determined by adjusting 
the calculated positions of the GT giant resonances for ^®Ca, ""Zr and ^°*Pb. This gives a mass dependence with 
fi = 0.7. The same mass dependence is assumed for the pp force and the coefficient is determined by a fitting procedure 
to /3-decay half-lives of nuclei with Z < 40. The result foimd in ref. [7] is x^t = ^.2 /A" "^ MeV and k^Jj, = 0.58 /A^-"^ 
MeV. A word of caution is in order concerning this parametrization of the residual forces. It serves to our purpose 
of comparing the effects of different deformed mean fields on the GT strength distributions, but one should keep in 
mind that the coupling strengths obtained in this way depend in particular, on the model used for single particle 
wave functions and on the set of experimental data considered. In ref. [7] a Nilsson potential was used and the set 
of experimental data did not include the nuclei under study here. Therefore, the coupling strengths of ref. [7] cannot 
be safely extrapolated and are not necessarily the best possible choices. As we shall see in the next sections, the 
strengths from ref. [7] reproduce well the data when using the WS potential, but one needs a somewhat smaller value 
of to reproduce the GT resonance with the HF mean field. 

The proton-neutron quasiparticle random phase approximation phonon operator for GT excitations in even-even 
nuclei is written as 

= E + , (2.9) 

TTI/ 

where a"*" (a) are quasiparticle creation (annihilation) operators, ojk are the RPA excitation energies, and X!^^ ,Y^jf 
the forward and backward amplitudes, respectively. The solution of the QRPA equations can be found solving first a 
dispersion relation [21], which is of fourth order in the excitation energies ujk- 

In the intrinsic frame the GT transition amplitudes connecting the QRPA ground state ]0) (P^^^ )0) =0) to one 
phonon states \cok) {^i^ |0) = \i^k)), are given by 

{uK\cTKt^\Q) =TM^'' . (2.10) 

The functions can be found for instance in [21]. The basic ingredients in their structure are the spin matrix 

elements connecting neutron and proton states with spin operators 
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^7 = {vWk\t^) . 

which can be written in terms of the coefRcients of the expansion in eqs.(2.1)-(2.2). 



K 



(2S) Vl + I^l 



(2.11) 

(2.12) 
(2.13) 



Once the intrinsic amphtudes are calculated according to eq. (2.10), the GT strength B{GT)± in the laboratory 
system for a transition /^/^^(O+O) — > IfKf{l~^Kf) can be obtained as 



2 

B{GT)± = ^ [Sk,,o {<i>K, |aoi±| 0o)' + 2<5k,4 {'t'Kj \o^t^\ 



(2.14) 



where wc have used the initial and final states in the laboratory frame expressed in terms of the intrinsic states \(I)k) 
using the Bohr-Mottelson factorization [22]. 

In the simple uncorrelated 2qp approximation, neglecting the residual ph and pp forces, the functions M^^ reduce 
to the following expressions 

= u^v^Y."^; M^'' = v^u^Y."^, (2.15) 

where the excitation energies are the bare two quasiparticle energies = E^, + E,^. 
The Ikeda sum rule is always fulfilled in our calculations 

Y^\{M'^f -{MXf'\='i{N -Z). (2.16) 



III. BULK PROPERTIES 

In this Section wc present results for the bulk properties of ^*'Gc and ^^Sc obtained from WS and HF descriptions. 

First, we analyze the energy surfaces as a function of deformation. In the case of WS, this is simply done by 
varying the quadrupole deformation of the potential /?2, which is an input parameter. In the case of HF, we perform 
constrained calculations [23], minimizing the HF energy under the constraint of keeping fixed the nuclear deformation. 

We can see in fig. 1 the total energy plotted versus the microscopically calculated mass quadrupole moment. The 
results correspond to HF calculations with the forces SG2 (solid line) and Sk3 (dashed line) , as well as to calculations 
with the WS potential (dotted line). The origin of the energy axis is different in each case but the distance between 
ticks corresponds always to 1 MeV. 

We observe that the HF calculation predicts the existence of two energy minima close in energy, giving rise to shape 
isomers in these nuclei, while the WS potential originates a single energy minimum, which is in agreement with the 
absolute prolate minimum in the case of ^^Ge and close to the prolate HF solution in the case of ''^Se. 

We can see in Table 1 the experimental and the microscopically calculated charge root mean square radii r,., 
quadrupole moments Qp, and quadrupole deformations [3 {(i — y^-§:p-)- In the case of ""^Se, the calculated values 
correspond to prolate/oblate deformations. The input WS prolate deformation is chosen to be (32 — 0.10 in both 
nuclei '^^Ge and ''^Se. In the oblate case of the nucleus ""^Se, the WS deformation chosen is (32 = —0.20. With these 
values we guarantee that the intrinsic deformations of the ground state are similar in HF and WS and therefore the 
differences in their predictions will have their origin in the structure of mean fields having the same deformation. 

The values obtained for the charge radii are in good agreement with the experimental values from rcf. [24], which 
are also shown in Table 1. They are also in good agreement with the results obtained from relativistic mean field 
calculations [25]: rc(rer)C^Ge) = 4.057 fm and rc(rei)('^'^Se) = 4.119 fm. 

The charge quadrupole moments quoted in Table 1 have been derived microscopically from the deformed potentials 
as ground state expectations of the Q20 operator. We can compare again with the results from relativistic mean field 
calculations of ref. [25]: Qp{Tei)C^Ge) = 111.4 fm^ and (5p(rei)(^^Se) = —146.8 fm^. These relativistic results are in 
perfect accordance with our calculated results. They can also be compared with experimental intrinsic quadrupole 
moments from ref. [26]. The empirical intrinsic moments are derived from the laboratory moments assuming a 
well defined deformation. These values are shown in Table 1 in the first place: Qp(exp)C^Ge) = 66(21) fm^ and 
Qp(cxp)C^Se) = 119(25) fm^. Experimental quadrupole moments can also be derived [27] from the experimental values 
of I3{E2) strengths, although in this case the sign cannot be extracted. Assuming that the intrinsic electric quadrupole 
moments are given by Q = ^/16^vB{E2)/5e^^, then \Qp(e^p)\C^Ge} = 164(24) fm^ and |Qp(exp)|(^^Se) = 205(24) fm^. 
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IV. GAMOW-TELLER STRENGTH DISTRIBUTIONS 



In this Section we show and discuss the Gamow- Teller strength distributions obtained from different choices of the 
deformed mean fields and residual interactions. 

As a general rule, the following figures showing the GT strength distributions are plotted versus the excitation 
energy of daughter nucleus. The distributions of the GT strength have been folded with Breit-Wigner functions of 
1 MeV width to facilitate the comparison among the various calculations, so that the original discrete spectrum is 
transformed into a continuous profile. These distributions are given in units of g\/4Tr and one should keep in mind 
that a quenching of the gA factor, typically gA,eff = (0.7 — 0.8) 5^, free, is expected on the basis of the observed 
quenching in charge exchange reactions. 

First of all, we discuss in figs. 2 and 3, the dependence of the GT strength distributions on the deformed quasiparticle 
mean field of ^^Ge and ^®Sc, respectively. To make the discussion meaningful we show the results obtained at the 
two-quasiparticle level without including the spin-isospin residual interactions. In these figures we can see the B{GT-) 
and B{GT+) strength distributions in the upper and lower panels, respectively. One should notice that the relevant 
strength distributions for the double /3-decay of ^^Ge, as it can be seen schematically in fig. 4, are the i?(GT_) 
distribution of the parent ^®Ge and the B{GT^) distribution of daughter ^^Se, but for completeness we show both 
distributions for each nucleus. Solid lines in figs. 2 and 3 correspond to the results obtained from the Skyrme force 
SG2 within a HF scheme, dashed lines correspond to the results obtained with the WS potential. The deformation of 
the mean fields are as indicated in Table 1, using the prolate shape in ''^Se. Pairing correlations are included in HF 
and WS cases in a similar way with the gap parameters for neutrons and protons mentioned earlier. Then, the only 
source of discrepancy between HF and WS comes from the different single particle wave functions and energies. 

In general, we observe that WS and HF produce a similar structure of three peaks in the B{GT ) profiles of ''^Ge 
and ■^'^Se, although the WS result s are somewhat displaced to lower energies with respect to the HF peaks. The 
strengths contained in the peaks are also comparable. In the case of the B{GT^) distributions, we first observe the 
different scale, which is about one order of magnitude lower than the B{GT-) scale. This is a consequence of the 
Pauli blocking. We can see from eq. (2.15) that while the occupation amplitudes u's and v' s favor M_ strengths, they 
are very small factors in M_|_ strengths when connecting similar proton and neutron states. The difference between 
total B(GT-) and B{GT+) strengths (Ikeda sum rule (2.16), which is fulfilled in our calculations), is a large number 
3(iV - Z) = 36 in ^^Ge and 3{N - Z) = 2A in ^^Se, reflecting the different magnitude of the B{GT_) and B{GT+) 
strengths shown in figs. 2 and 3. 

The profiles of the B{GT+) distributions with WS and HF present some discrepancies that are amplified by the 
scale. In particular, it is remarkable the large strength produced by WS in the region of high excitation energies in 
''^Ge that we discuss later in terms of the single particle wave functions. 

In order to clarify the origin of the various peaks in the strength distributions we have added in fig. 2 labels showing 
some of the leading transitions generating the strength. The labels stand for pK^ — nK'^ of the orbitals connected 
by the spin operator in eq. (2.11) and a number that identifies the transition. In both cases, B{GT-) and B{GT+), 
the same type of transitions are connected by the GT operator but the occupation probabilities, weighting the matrix 
elements, enhance or reduce them accordingly. We can see from fig. 2 that the structure of the profiles in both WS 
and HF are generated by the same type of GT transitions. 

This can be further illustrated by looking at fig. 5, where we show the single particle energies for protons and 
neutrons obtained in HF(SG2) and WS in ^®Ge. In the left part of the figure corresponding to the HF calculation we 
have plotted the occupation probabilities and and the Fermi energies X„ and Att- We can also see for completeness 
the spherical levels labeled by their ij values. We have indicated by arrows the most relevant Gamow- Teller transitions 
in the f3~ and /?+ directions that are labeled by the same numbers used in fig. 2 to identify the peaks. To be more 
precise, we can see in Table 2 the correspondence between these labels and the transitions connecting the proton and 
neutron states using the asymptotic quantum number notation [A^n^AJif . 

Now, looking at fig. 2, we can understand that the two first peaks in B{GT-) are generated mainly by transitions 
between neutrons and protons dominated by contributions within the N = 3 shell and that the third peak is generated 
by transitions between neutrons and protons with main contributions coming from the iV = 4 shell. The different 
energies of the peaks are due to the different concentration of energy levels in HF and WS. 

In the case of B{GT+), the strength below 8 MeV is mainly generated by transitions within the N = 3 shells. 
Beyond 8 MeV the strength, which is negligible in HF, is generated by transitions between the proton shell N = 2 and 
the neutron shell iV = 4 as well as between the proton shell N = 3 and the neutron shell N = 5, always understood 
as the main components of the wave functions. Then, very deep inside protons {Vp = 1) are connected with very 
unoccupied neutron states (u„ = 1), giving rise to maximum occupation factors. The different behavior in this high 
energy region between HF and WS is therefore due, other factors like deformation and occupations being equal, to 
the structure of the deformed orbitals. 
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To illustrate the role of the different single particle wave functions in the development of the peak structure, we 
consider in detail the case of the last peak observed in the B{GT^) distribution of the WS potential. As we can see 
it is mainly due to a transition between the proton state [303] in the = 3 shell with K"^ = 7/2~ and the neutron 
state [523] in the N = 5 shell with K'^ = 5/2~. The structure of the single particle wave functions, according to the 
expansion in eq. (2.1), of these two states in the cases of HF and WS can be seen in Table 3. With these coefficients 
we can construct the matrix elements in eq. (2.12). The resulting strength is almost two orders of magnitude in favor 
of WS, which explains the huge discrepancy observed between WS and HF in the higher energy domain. 

Nevertheless, these discrepancies are smaller in the case of the B{GT^) of '^^Se, which is the relevant branch for 
the double /3-decay of the parent nucleus ^^Ge. 

Figures 6 and 7 contain the strength distributions obtained from QRPA calculations for ''^Ge and ''^Se, respectively. 
The data in Fig. 6 are from ref. [28] and were obtained from charge exchange ^^Ge(p,n)^^As reactions. The thick line 
in Fig. 6 corresponds to these data folded by the same procedure used for the theoretical results. The data in Fig. 7 
and 8 are from ref. [29] and were obtained from charge exchange ^^Se(n,p)''^As reactions. 

The coupling constants of the ph and pp residual interactions used in Figs 6-8 are from ref. [7] in the case of WS. In 
our case with A = 76, these parameters are Xg-t ~ '-'•25 MeV and K^fj, = 0.027 MeV. In the case of the HF calculations 
with the Skyrme forces Sk3 and SG2 better agreement with the measured location of the B{GT-) resonance in ^^Ge 
is obtained with a somewhat smaller value of the ph strength. The curves shown in Figs. 6-8 for the HF results have 
been obtained using x^'jp = 0.16 MeV and the same kF^j, = 0.027 MeV. 

In Fig. 6 we have used the prolate deformations for ""^Ge given in Table 1. We can see that WS follows the 
structure of the experimental B{GT-) strength distribution with two peaks at low energies {Eex =5 and 8 MeV) and 
the resonance at 11 MeV. The HF calculations produce also a few peaks at low excitation energies and a resonance 
between 10 and 13 MeV. We can see that the structure of the strength distributions is qualitatively similar for the 
two Skyrme forces and that the difference with the WS curves can be traced back to the discrepancies found at the 
two quasiparticle level. 

Fig. 7 contains similar calculations for ^^Se. The coupling strengths of the residual forces are as indicated for ^^Ge. 

The results in the HF cases are obtained with the oblate deformation of ^''Sc that produces the absolute minimum of 
the energy and agrees better with the experimental quadrupole moment. In general, comparison with experiment is 
reasonable and should not be stressed too much since, as stated in ref. [29], the experimental results, especially above 
6 MeV, must be considered to be of a qualitative nature only. 

The role of the residual interactions on the GT strengths was already studied in ref. [8], where it was shown that 
the repulsive ph force introduces two types of effects: A shift of the GT strength to higher excitation energies with 
the corresponding displacement of the position of the GT resonance and a reduction of the total GT strength. The 
residual pp, being an attractive force, shifts the strength to lower excitation energies, reducing the total GT strength 
as well. Also shown in ref. [8] was the effect of the BCS correlations on the GT strength distribution. The main effect 
of pairing correlations is to create new transitions that are forbidden in the absence of such correlations. The main 
effect of increasing the Fermi diffuseness is to smooth out the profile of the GT strength distribution, increasing the 
strength at high energies and decreasing the strength at low energies. 

The role of deformation was also studied in ref. [8], showing that the GT strength distributions corresponding 
to deformed nuclear shapes are much more fragmented than the corresponding spherical ones, as it is clear because 
deformation breaks down the degeneracy of the spherical shells. It was also shown that the crossing of deformed energy 
levels that depends on the magnitude of the quadrupole deformation as well as on the oblate or prolate character, 
may lead to sizable differences between the GT strength distributions corresponding to different shapes. 

We can see in fig. 8 the GT strength distributions in '^^Se obtained from spherical, prolate and oblate shapes. They 
correspond to QRPA calculations performed with the HF basis obtained with the force Sk3. In the spherical case, 
the only possible transitions (see fig. 5) are those connecting spherical partners with Ai? = 0, A j = 0, 1, in allowed 
approximation. Therefore, there is GT strength only at a few excitation energies. The strength we observe in fig. 8 
is the result of the folding procedure performed at these energies. On the other hand, in the deformed cases we can 
observe a stronger fragmentation, which is the result of all possible connections among the deformed states (see fig. 
5). Thus, the spherical peaks become broader when deformation is present. 

We can see in Table 4 the total GT strengths in ^^Se contained below an energy cut of 60 MeV. We show the results 
obtained for /?+ and j3~ strengths with oblate, spherical, and prolate shapes. The Ikeda sum rule 3(A'' — Z) = 24, is 
fulfilled at this energy cut within a 0.3% accuracy. We can see from Tabic 4 that deformation increases both /?+ and 
/3~ strengths in a similar amount in order to preserve the Ikeda sum rule (/3^ — /?+). We also show for comparison the 
results obtained in 2qp approximation. We can see the reduction of the strength introduced by the QRPA correlations, 
which is again similar in absolute terms for /3+ and strengths in order to keep the Ikeda sum rule conserved in 
QRPA. Since the fi~ strength is much larger than the /3+ strength, the relative effect of the QRPA correlations is 
much stronger for /?+, where the total strength is reduced by a 50%. 
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Comparing the results for ^^Se obtained at different deformations with the selfconsistent mean fields (HF with 
Sk3) in Fig. 8 and Table 4, we see that there is a strong dependence on deformation in the strength distributions 
as a function of the energy. However, the total strength does not depend so much on deformation. There is an 
increase of a few percent in going from the spherical to the oblate and prolate shapes. The latter observation enters in 
contradiction with SU(3) and shell model calculations by previous authors [30] on the dependence on deformation of 
the GT strengths in ^"Ne and ^^Ti. We think that this is due to the much larger and richer single particle basis used 
in the present calculations. In our case each single particle state contains mixtures from many harmonic oscillator 
shells (up to N = 10), while in the above mentioned calculations [30], the single particle basis is restricted to a 
single harmonic oscillator major shell (the sd shell in ^*^Ne and the fp shell in ^^Ti). On the other hand, one may 
question whether in the deformed cases the total strengths calculated here may contain spurious contributions from 
higher angular momentum components in the initial and final nuclear wave functions. Since the matrix elements of 
the transition operator, which is a dipole tensor operator, are taken between the states considered in the laboratory 
frame (see eq. (2.14)), the effect of angular momentum projection is to a large extent taken into account. We have 
calculated an upper bound to such contributions using angular momentum projection techniques [31]. We find that 

this upper bound is less than 1% (^ (j^) , with {J]_) = 19 for the oblate shape in ^''Sc). Thus, exact angular 
momentum projection would not wash out the small increase of the total strength with deformation. 

V. CONCLUDING REMARKS 

We have studied the GT strength distributions for the two decay branches f3~ and /?+ in the double /9-decay of ""^Ge. 
This has been done within a deformed QRPA formalism, which includes ph and pp separable residual interactions. 
The quasiparticlc mean field includes pairing correlations in BCS approximation and it is generated by two different 
methods, a deformed HF approach with Skyrme interactions and a phenomenological deformed WS potential. One 
difference is that with HF and Sk3 we get the minimum and stable deformation for ''^Se to be oblate, while the prolate 
minimum is comparable to that obtained with WS and it is higher in energy. 

We have studied the similarities and differences observed in the GT strength distributions with these two methods. 
Among the similarities we can mention the structure of peaks found in the strength distributions and among the 
differences the displacement in the excitation energies found between HF and WS results. This discrepancy has its 
origin in the structure of the single particle wave functions and energies generated by the deformed mean fields. 
This also implies that different mean fields require different residual interactions to reproduce the experimental GT 
resonances. 

Therefore, in order to obtain reliable GT strength distributions and consequently reliable estimates for double 
/3-decay half-lives, it is important to have not only the proper residual interactions but also a good deformed single 
particle basis as a starting point. In the case of HF we have seen that standard Skyrme forces, such as SG2 or 
Sk3, give a good description of the GT strength distributions provided the proper residual interactions are included. 
Even though the selfconsistent HF approach is a more sophisticated type of calculation, the deformed WS potential 
produces comparable results when the parameters of the potential and the residual interactions for a given mass region 
are chosen properly. 

There is work in progress to extend these calculations to the double /?-decay process studying the dependence on 
deformation of the half-lives. 
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Table 1. Charge root mean square radii Tc [fm], intrinsic charge quadrupole moments Qp [fm^], and quadrupole 
deformations j3 for ^®Ge and ^®Se calculated with various assumptions for the mean field. In the case of ''^Se we 
show theoretical values corresponding to the prolate shape in first place and to the oblate shape in second place. 
Experimental values for Tc are from [24] and for Qp from [26] the first value and from [27] the second (see text). 









Qp 


P 




exp. 


4.080 - 4.127 


66(21) - 164(24) 


0.10 - 0.24 




Sk3 

SG2 
WS 


4.130 

4.083 

:\sm 


111.0 

105.9 
110.9 


0.161 
0.157 

().17() 




exp. 


4.088 - 4.162 


119(25) - 205(24) 


0.16 - 0.29 




SkS 
SG2 
WS 


prol / obi 
4.170 / 4.180 

4.113 / 4.143 
3.991 / 4.138 


prol / obi 
117.5 / -136.0 
35.2 / -140.6 
81.6 / -141.4 


prol / obi 
0.158 / -0.181 
0.049 / -0.191 
0.119 / -0.193 



Table 2. Correspondence of the labels used in figs. 2 and 5 with the asymptotic quantum numbers notation 





P- 


/3+ 


(1) 


:/[301]l/2- ^ 7r[301]3/2- 


TT [303] 7/2" ^ 1/ [303] 5/2" 


(2) 


i/[301]3/2- ^ 7r[301]l/2- 


7r[312]5/2- ^ !y[303]5/2- 


(3) 


i/[303]7/2~ ^ 7r[303]5/2- 


7r[312]5/2- !/[312]3/2- 


(4) 


!^[312]5/2- ^ 7r[312]3/2- 


7r[202]3/2+ ^ j/[413]5/2+ 


(5) 


[420] 1/2+ TT [440] 1/2+ 


7r[330]l/2- i^[530]l/2- 


(6) 




7r[303]7/2- ^ i^[523]5/2- 



Table 3. Main coefficients in the expansion of eq. (2.1) for the proton state [303] with K'^ = 7/2^ and the 
neutron state [523] with = 5/2~. This is the main contribution to the peak at 15 MeV in the i?(GT+) strength 
distribution of ^^Ge with the Woods-Saxon potential. The basis states are labeled by | A^^n^ A) quantum numbers. The 
table contains also the contributions from these basis states to the spin matrix elements in eqs. (2.11)-(2.12). 







|303) 


|503) 


|523) 


|703> 


|723) 


|903) 


7/2~ proton 


















HF(SG2) 


-0.9742 


0.2204 


-0.0061 


0.0219 


-0.0272 


0.0122 




WS 


0.9876 


-0.1400 


0.0563 


-0.0233 


0.0107 


-0.0295 


5/2~ neutron 


















HF(SG2) 


0.1369 


0.5933 


-0.5031 


-0.3928 


0.2349 


0.2385 




WS 


-0.2397 


-0.5173 


0.5049 


0.3794 


-0.2596 


-0.2056 


contribution to 


















HF(SG2) 


-0.1333 


0.1308 


0.0031 


-0.0107 


-0.0064 


0.0029 




WS 


-0.2367 


0.0724 


0.0284 


-0.0088 


-0.0028 


0.0061 
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Table 4. Total Gamow- Teller strength in ^®Sc calculated with the force Sk3. Results correspond to /?+ and /3~ 
strengths for the oblate, spherical, and prolate shapes calculated in 2qp and QRPA approximations. All the GT 
strength contained below an excitation energy of 60 MeV has been included. 







oblate 




spherical 




prolate 




/3+ 


0- 


/3+ 


P- 


/3+ 


P- 


RPA 


2.420 


26.331 


1.846 


25.765 


2.599 


26.524 


2cip 


4.387 


28.298 


3.816 


27.736 


4.971 


28.892 
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FIG. 1. Total energy as a function of the mass quadrupole moment obtained from deformed Hartree-Fock calculations with 
the Skyrme forces SG2 (solid line) and Sk3 (dashed lines), and from deformed Woods-Saxon potentials (dotted line). The origin 
of the energy axis is different in each case but the distance between ticks corresponds always to 1 MeV. 
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FIG. 2. Gamow- Teller B{GT^) and B{GT+) strength distributions [.g_4/47r] in '^''Gc plotted as a function of the excitation 
energy of the daughter nucleus. We compare results of HF(SG2)+BCS (solid line) and WS+BCS (dashed line) approximations 
for the prolate minima. 
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FIG. 3. Same as in fig. 2 for ^«Se. 
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FIG. 6. Gamow-Teller B(GT-) and B{GT+) strength distributions [g'i/4Tv] in ™ Gc plotted as a function of the excitation 
energy of the daughter nucleus. We compare QRPA results of HF(SG2) (solid lines), HF(Sk3) (dotted lines) and WS (dashed 
lines). Deformations and coupling strengths of the residual interactions are given in the text. Experimental data (thick solid 
lines) are from [28]. 
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FIG. 7. Same as in fig. 6 for ''^Se. Data are from [29]. 
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FIG. 8. QRPA Gamow- Toller B(GT-) and B(GT+) strength distributions in '^^Se. The ealculations correspond to the force 
Sk3 for spherical (dotted line), prolate (solid line) and oblate (dashed line) shapes. Data are from [29]. 
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